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Quantum Opacity, the RHIC HBT Puzzle, and the Chiral Phase Transition 
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We present a relativistic quantum mechanical treatment of opacity and refractive effects that 
allows reproduction of observables measured in two-pion (HBT) interferometry and pion spectra at 
RHIC. The inferred emission duration is substantial. The results are consistent with the emission 
of pions from a system that has a restored chiral symmetry. 

PACS numbers: 25.75.-q 
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The space-time structure of the "fireball" produced in 
the collision between two heavy ions moving relativis- 
tically is studied by measuring the two-particle momen- 
tum correlations between pairs of identical particles. The 
quantum statistical effects of symmetrization cause an 
enhancement of the two-boson coincidence rate at small 
momentum differences that can be related to the space- 
time size of the particle source. This method, called HBT 
interferometry, has been applied extensively in recent ex- 
periments at the Relativistic Heavy Ion Collider (RHIC) 
by the STAR and PHENIX collaborations Q- 

The invariant ratio of the cross section for the pro- 
duction of two pions of momenta Pi,P2 to the prod- 
uct of single particle production cross sections is ana- 
lyzed as the correlation function C(pi,p2). We define 
q=Pi~P2 and K=(pi+p2)/2, with as the compo- 
nent perpendicular to the beam direction. (We focus 
on mid-rapidity data, where K = Kt.) The correla- 
tion function is parameterized for small q as C(q, K) 
«1 + A(l — R 2 D q 2 - R 2 s q 2 s — Rfqf),where o, s, I represent 
directions parallel to Kt, perpendicular to Ky and the 
beam direction, and parallel to the beam direction^. 
Early Q and recent Q hydrodynamic calculations pre- 
dicted that a fireball evolving through a quark-gluon- 
hadronic phase transitions would emit pions over a long 
time period, causing a large ratio R /R s . The puzzling 
experimental result that R Q /R S m 1 ispart of what 
has been called "the RHIC HBT puzzle" Q. 

Data shows the medium produced by 200 GeV Au+Au 
collisions to be very dense. Consequently, pions should 
emerge from an opaque source 0. Analyticity tells 
us that opacity implies accompanying refractive effects. 
Our purpose here is to derive and apply a relativistic 
quantum-mechanical treatment of opacity and refrac- 
tive effects that simultaneously reproduces the values 
of R a ,R s ,Ri and the pion spectrum for central RHIC 
Au+Au collisions I 



ization with approximately boost-invariant longitudinal 
dynamics: 



The experimental observables depend on an emission 
function, the Wigner transform of the density matrix for 
the currents that emit pions. This emission function has 
often been modelled (see review |j|), as a function So 
having the form of a hydrodynamic source parameter- 
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in which components of x^ are expressed using the vari- 
ables t = \Jt 2 — z 2 , ?y = ^logf^f, b = (xi,X2), and 
K ■ u = Mt cosh 77 cosh r\ t (b) — Kt sinh r\ t (b) cos 4>, where 
4> is the angle between Kt and b, Mt = \/ K' T + m 2 , 
and [i-k is the pion chemical potential. The function p(b) 
represents the cylindrically symmetric transverse source 
density. We use p{b) = [l/(exp((r — Rws)/ a ws) + I)] 2 , 
reflecting superimposed nuclear densities. Eq. rep- 
resents the transverse flow rapidity using a linear radial 
profile of strength r\f. rj t {b) = r]f R . Eq. Q incorpo- 
rates the finite lifetime and size of the source: pions are 
emitted for a duration controlled by the parameters At 
and A?y. Using Eq. Q is not sufficient to capture all of 
the physics. As the basis of the blast wave parameteriza- 
tion, it gives At f=a and does not predict the magnitude 
of the pion spectrum |lf|. 

The salient feature of the 200 GeV data is the high den- 
sity of the produced matter, so we treat the effects of pion 
interactions with the dense medium. We adopt a single- 
channel approach that uses the interaction-distorted in- 
coming wave ^^{xi) in which[ll|: 

S(x,K) = J d±K>S (x,K>) J 0^e-^'-*' (4) 

x ^l-\x + x'/2)^-J*(x-x'/2). 

One could imagine a generalization in which would 
have two or more components (corresponding e.g. to a 
p final state in addition to a pion final state). Our task 
here is to compute pion correlation functions, so that it 
is sufficient to use a single channel wave function. One 
obtains the single-pion emission function from Eq. Q by 
using the same momentum (either pi or P2) to compute 



the wave function ^ 
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Using Eq. requires evaluating an eight-dimensional 
integral, and modelling the interactions that determine 
^pi • We use symmetries to reduce the number of inte- 
grals and obtain a tractable treatment of the interactions. 
First, note that (re) is an energy-eigenfunction [ill ]: 

*p (a;) = e~ iUp x °^^'(-x). We assume that the matter 
formed in the central region of the collision is cylindri- 
cally symmetric with a very long axis, so that 

*(-> 2 (x) = e=F'«.*/V(-)(x J . = b) > (5) 
pi.2 = K±q T /2±2q i /2, 

with "i/ip ■'(b) obtained by solving a two-dimensional 
Klein-Gordon equation 

(-Vi +U(b)) (b) =p%->(b). (6) 

The "optical potential" U is a complex, azimuthally- 
symmctric function depending on pion momentum and 
local density that represents the strength of the inter- 
action between a pion and the medium. Within our for- 
malism the influence of time-dependent effects in U intro- 
duced by the time-dependent source So is incorporated 
in the energy dependence of the optical potential. How- 
ever, we note that the pions-medium interaction time is 
restricted by So. 

At large values of K the solution of Eq. © reduces to 
well-known semi-classical (eikonal) expression, but our 
procedure is more general. For example at K=0 there 
is no distinction between the out and side directions, so 
Ro=R s - This constraint is violated if the eikonal approx- 
imation to © is used. We use a partial wave expansion: 

4">*(b) = f {p, 6) + 2 fm(P, b)(-i) m cosmcp (7) 

m— l.oc 

to solve Eq.(0 exactly and maintain the constraint. 

The optical potential accounts for situations in which 
the pion changes energy or disappears entirely due to its 
interactions with the dense medium. Suppose, e.g., that 
the medium is a gas of pions. Then irir scattering would 
be the origin of U . In the impulse approximation, the 
central optical potential would be Uq = —Airfpo, where 
/ is the complex forward scattering amplitude and po 
the central density. For low energy pion-pion interac- 
tions, 47rlm[/(p)] = pa, with a w 1 mb. At a mo- 
mentum p = 1 fin - , using a pion density about ten 
times the baryon density of ordinary nuclear matter, 
Im[E/(0)] ks —0.15 fm~ , representing significant opac- 
ity. Furthermore, if two interacting pions each have less 
energy than half of the rho meson mass, the final state 
interactions would be strongly attractive. 

We model the interaction U for a situation in which the 
medium is dense, but we have not otherwise specified its 
nature. Now we consider chiral symmetry, which gives 
a general form for the dispersion relation of low energy 



pions in nuclear matter: 

co 2 = v 2 (f+ml(T)). (8) 

Here m ff (T) is the ir screening mass, v is the so-called tt 
velocity, p is the momentum operator, and the pole mass 
is m OT (T)u. Son and Ster)hanov|l2j argued that v and the 
pion pole mass decreases at temperatures near the critical 
temperature, while the screening mass increases. Several 
papers have used specific low-order model calculations 
that challenge this result fl3l H3| . Here we start with the 
general form |JSJ that appears in both Refs. and 
use data to determine the parameters. We first define an 
equivalent U by using J 1 = p 2 + m 2 , Eq. (JSJ, and the 
Klein Gordon equation p 2 + U = p 2 , to obtain: 

U=-{ml-v 2 ml{T))-{l-v 2 )p 2 . (9) 

For matter of finite size, Eq. @ suggests a potential 
with terms constant and proportional to p 2 . The latter 
has the form W2P 2 — V ■ w^V, where w 2 .k are complex 
constants. For the 200 GeV data, the results prefer a 
small value of wx, so we take wk = 0. The optical 
potential is therefore modelled as: 

U p {b) = -{w + w 2 p 2 )p{b), (10) 

with wq real (no opacity at p=0). The density profile 
p(b) is that specified in Sq(x, K). This simple form is 
sufficient to account for the data we study. 

Eans. (|lllQ|) are the essence of our model. These may 
be compared with the Buda-Lund model, which is an 
efficient representation of the data|l5j in which the tem- 
perature and fugacity are taken as position-dependent 
functions appearing in a Boltzmann distribution. The 
effects of our optical potential could provide an explana- 
tion of those deduced dependencies. 

The use of Eqs. loll I) in Eq. (@J gives 

S{x,K) = -I-SoOr,^ *-*** J d 2 b'B n (b,b') 

x ^ p ;>(b + b'/2)<A p ;>*(b-b72), (ii) 

where B,(b,b') = / d 2 K' T B v (b, K' T ) exp [—iK' T ■ b'} . 
The range of the variable b' appearing in ^(bjb') is 
controlled by a size, 1/T that is much smaller than the 
source size. Thus it is reasonable to ignore the ±b'/2 
appearing in (|ll|l . Maintaining the phases is important, 
so we use V P _) *(b±b72) = e- i P-( b±b '/ 2 ) Xp (b±b'/2) -> 
e -tp-(b±b / 2 )^ p (b). This approximation is exact in sev- 
eral different limits: plane wave, short wavelength, long 
wave length, and taking (RwsT — > oo). 

We expand the exponentials in l|ll(l involving q°, qi to 
2nd order, expand the Bose-Einstein function in a series 
of Boltzmann functions of temperature T n = T/n, and 
analytically evaluate integrals over r and 77 to obtain: 

C(q, K T ) = 1 - q ?R? q l0 2 Ar 2 + (12) 
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where £„(6) = M T cosh Vt {b)/T n + 1/A V 2 , /i(0 = 
2K (0/£ + 4Ki(0/£ 2 , /o(0 =■ 2^(0, / 2 (0 = ^o(0 + 
#a(0> /s(0 = 2(^(0 + ^8(0/0. and ^ are modified 
Bessel functions. The pion spectrum is given by: 
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The angular integrals are performed analytically 
E q. 0. The tr ansverse HBT radii are R t {K T ) = 
yjl - C(Aqi,K T )/ Aq i7 with i = o,s and Aqi « if T /40. 

The values of the parameters T,rjf, At, RwSidws^ 
100,2, To, At? and /i w arc varied to reproduce the STAR 
200 GeV data for R , R s , Ri Q and the magnitude and 
shape of the pion spectrum [8j . The resulting parameters 
are displayed in Table I, with the variances that produce 
a x 2 increase of one unit. The agreement of the curves 
with STAR HBT radii Q in Fig. 1 and with the STAR 
pion spectrum || in Fig. 2 is quite good, giving a. x 2 ~ 
3.7 per data point and 5.6 per degree of freedom. 

We need to assess and interpret the values of our pa- 
rameters. The temperature (T— 173 McV) is comparable 
to the T c expected for a chiral phase transition (wl60 
McV). The flow rapidity (77/= 1 .31) corresponds to a max- 
imum flow velocity of 0.85 c, which is not unreasonable. 
The source size is large (i?ivs=H-7 fm), 4.4 fin larger 
than a gold nucleus (R=7.3 fm). If the system requires an 
expansion time tq= 8.22 fm/c to expand by this amount, 
the average expansion velocity would be about 0.5 c, 
which is reasonable. The emission duration (Ar=2.9 
fm/c) is comfortably smaller than To- The longitudinal 
width (Ar]= 1.063) is smaller than expected, but is large 
enough that the system's axial length (~ 2tqAti ~ 17.5 
fm) is sufficient for the approximate validity of Eq. (J5J . 

We assess the strength of U by examining it at p w 1 
fin . At the system center, the momentum-dependent 
term of U is -(0.58+0.12i) fm~ 2 . To understand the 
imaginary part, which is comparable to our estimate 
given above (and consistent with the presence of high 
density matter), consider the equivalent classical mean 
free path: p/|/m[£/]| « 8 fm <C 2Rws- Thus, the imagi- 
nary potential is large enough to restrict emission of pi- 
ons from regions deep inside the medium. To understand 
the real part, note that the strength of the attraction is 
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FIG. 1: (Color online) HBT Radii R a ,Ro„ Ri and the ratio 
Ro/R s \ Data 0): V (green) => A (red) =>■ -k'tt~. 

Curves: solid (black) =>■ full calculation; dotted (green) 
r/f — (no flow); dashed (red) => Re[(7]=0 (no refraction); 
dot-dashed (blue) =>■ (7=0 (no potential), double-dot-dashed 
(magenta) =>• substituting Boltzmann for Bose-Einstein ther- 
mal distribution. Insets show predictions of \ow-Kt resonance 
behavior in R and R 3 . 



greater than m 2 = 0.49 fm -2 . Thus inside the medium, 
the pion acts as if it has no mass. This is why we assert 
that a chiral phase transition has occured. 

We further recall Eq. JjJJl, and note that a momentum- 
dependent term of -0.58 fm -2 corresponds v = 0.65. The 
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TABLE I: Parameters of the calculation with variances 



T(MeV) 


Vf 


Ar(/m/e) 


Rws(fm) 


aws(fm) 


w (fm z ) 


W2 


r (/m/c) 


A77 




173.2 


1.314 


2.852 


11.728 


0.725 


0.137 


0.582 ± i 0.121 


8.23 


1.063 


123.2 


±1.6 


±0.025 


±0.067 


±0.056 


±0.015 


±0.046 


±0.014 ±0.002 


±0.10 


±0.032 


±1.1 
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FIG. 2: 
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(Color online) Pion momentum spectrum. Data y|: 
) => 7r + ; A (red) => 7r _ . Inset is Iow-T^t prediction. 



momentum independent term of U (-.14 fm -2 ) corre- 
sponds to of vm n (T) = 0.6 fm -1 . These values are com- 
parable to the estimates of [13, but not too far from the 
lowest value v = 0.83 of Ref~[lJ|. 

The chemical potential and Bosc-Einstein distribution 
have modest effects on the radii but are very important 
for the normalization and slope of the pion spectrum. 

Figs. ^ |2 show the importance of both the real and 
imaginary parts of the optical potential: omitting cither 
drastically changes the predictions. The attractive real 
potential is the critical element needed to reproduce the 
Kt dependence of R and Rs- 

We predict peaks in R Q and R s at low momentum (p w 
15-65 MeV/c) and a rapid rise and peaking in the low- 
momentum spectrum. This is a pionic version of the 
Ramsaucr effect, in which cross sections show peaks when 
the scattered wave is in phase with the incident plane 
wave. We confirm the computed existence of such peaks 
through analytic calculations for a purely real attractive 
square well potential. The PHOBOS detector (and per- 
haps the BRAHMS detector) at RHIC could, in principle, 
confront these low-momentum predictions of structure. 

Our rclativistic quantum-mechanical treatment of re- 
fractive and opacity effects on two-pion correlations and 
7r spectra has surmounted a number of technical problems 
associated with previous semi-classical approximations, 
while achieving an excellent description of the 200 GeV 
STAR data. The results are consistent with the inter- 



pretation that the dense medium has undergone a chiral 
phase transition. The definitive predictions of interest- 
ing momentum dependences at small momenta should be 
testable in present and future RHIC experiments. 
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